Use of Lambert's Theorem for the n-Dimensional Coulomb Problem 
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We present the analytical solution in closed form for the semiclassical limit of the quantum 
mechanical Coulomb Green function in position space in n dimensions. We utilize a projection 
method which has its roots in Lambert's theorem and which allows us to treat the system as an 
essentially one dimensional problem. The semiclassical result assumes a simple analytical form and 
is well suited for a numerical evaluation. The method can also be extended to classically forbidden 
space regions. Already for moderately large principal quantum numbers v > 5, the semiclassical 
Green function is found to be an excellent approximation to the quantum mechanical Green function. 

PACS numbers: 03.65.Sq,03.65.Ge,03.65.Nk 



I. INTRODUCTION 

The laws of planetary motion remained for a long time 
a mind-puzzling challenge. It was Johannes Kepler who 
published 400 years ago his book Astronomia Nova which 
contained his famous first two laws on planetary motion. 
Kepler's conclusion that all planets move in elliptical or- 
bits with the Sun in one focus was based on his ingenious 
evaluation of very accurate observations of the path of 
the planet Mars by the astronomer Tycho Brahe, the 
last of many important astronomers who made their ob- 
servations without the help of a telescope. As is well 
known, the mathematical construction scheme for ex- 
pressing the motion of bodies in a gravitational (l/r) 
potential in mathematical terms goes back to the days of 
Newton's Principia Mathematica, first published in 1687. 
This work unifies Galileo's ideas about motion in a grav- 
itational field and Kepler's laws on planetary motion. 
In the 18*'' century it was still a major problem to follow 
the motion of a planet along its elliptical path or, more 
general, along a curved trajectory. For the 1/r-potential 
this difficulty was solved by the Swiss Mathematician and 
Physicist Johann Heinrich Lambert who proved geomet- 
rically that the transfer time along a planetary orbit con- 
necting two position vectors r and r' depends only on the 
two combinations a+ and 



and a_ 



(1) 



where s is the distance between r and r'. The position 
vectors are meant relative to the force center (in Lam- 
bert's case the Sun). The additional dependence of the 
travel time on E will be discussed later. Equation ([T]) 
is a peculiarity of the l/r-potential. The fact that the 
transfer time depends only on a+ and a_ is called Lam- 
bert's theorem. 



The agreement between the calculated and observed po- 
sitions of the planets was historically the most impor- 
tant success of classical physics. With the advent of 
quantum mechanics, the Kepler problem was replaced 
by the Coulomb problem for the hydrogen atom. Feyn- 
man's path integral method revealed the close connection 
between classical and quantum mechanics. The fixed- 
energy propagator for the Coulomb problem is known 
analytically both in configuration and momentum space 
0, 0, S, II i, H, S i, M- However the corresponding 
semiclassical approximation has not been given before in 
closed analytic form because of the appearance of a rather 
complicated prefactor, the so called Van Vleck-Pauli- 
Morette determinant. The semiclassical approach to the 
Coulomb problem in n > 1 dimensions the determinant 
has been calculated so far only numerically [l3|- Based 
on Lambert's theorem, we will derive a simple and use- 
ful analytic expression for the Van Vleck-Pauli-Morette 
determinant in n spatial dimensions. The result will put 
us in the position to derive a two-line expression for the 
semiclassical Green function. 



II. LAMBERT'S THEOREM FOR THE 
REDUCED ACTION 

It is a simple exercise in classical mechanics to analyze 
the relative motion for the Kepler or Coulomb Hamilto- 
nian 

-2 Kr 
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(2) 



where /z is the reduced mass and Kc the strength of the 
attractive l/r potential. The corresponding motion in 
time is given by 




dr . (3) 
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Here A — ^r'^ip is the angular momentum about the 
center of force for elliptic motion with semimajor axis 



a = Kc/{2\E\), for E<0. 

An important element for the transition from classical 
mechanics to quantum mechanics is the reduced action 
W, also called action integral S. In order to avoid confu- 
sion, we reserve S here for Hamilton's principal function 
(see below). Within the time-independent Hamilton- 
Jacobi theory the reduced action W is given by 



W{r,r';E) = p(f) • df 



(4) 



For elliptic motion in the x-y-plane an explicit evalua- 
tion of Eq. ^ is easily achieved by introducing for exam- 
ple Cartesian coordinates with the origin at the center of 
the ellipse 



a cos ^ , y = b sin <^ , 



(5) 



where b = a\/l — is the semiminor axis of the Kepler 
ellipse with eccentricity e. If we substitute Eq. (O into 
Eq. ([U and use 



i-t'-l/^(e-"ine-r + "inO 



(6) 



for the transfer time between two points on the ellipse, 
we obtain 



Wir,r';E) = y/^iaK^ (C + esin^ - - esin^')- (7) 

W{r,r';E) is a function of E and of the initial and fi- 
nal coordinates r and r' of the planet. Therefore other 
dynamical quantities, like the orbital angular momen- 
tum A must be eliminated. Hence we have to get rid 

0. 



of e = ^1 - 2\E\A^/{nK^) in Eq. 
manipulations (see App. |^ lead to 



A few algebraic 



W(r,r';E) = ^JfiaKc (7 + sin 7 - (5 - sin (5) (8) 



with 



o 7 r + r' + s 

r - ^ 

2 4a 



and 



4a 



(9) 



In the last equation and in what follows r and r' are the 
distances from the focus of the ellipse (i. e. the center 
of force) to two arbitrary points on the elliptical orbit. 
As before s stands for the distance between r and r'. 
The situation is depicted in Fig. [1] It was Lambert [llj 
who succeeded to map jdliptical motion to coUinear mo- 
tion. He also proved ([11], Lemma 24) that for fixed en- 
ergy E < Q, the classically allowed elliptic motion from 
a given initial point to a final point M can generally 
occur on two different ellipses unless we have circular mo- 
tion (compare [H, p. 27). From the last two equations it 
becomes obvious that the reduced action is a function 
of a+ = r + r' + s and a_ = r + r' — s. Note that the 
energy dependence of W enters through the semimajor 
axis a = Kc/{2\E\). 

Another piece of information is Hamilton's principal 




FIG. 1: Lambert's projection of elliptic motion to collinear 
motion. Shown are two ellipses with the same lengths of the 
semimajor axes ^AiBi = ^^2-62 and a common focus located 
at F. The centers of the two ellipses are denoted by Ci and 
C2. Lambert's lemma 24 allows to relate the motion from A'^ 
to M on both ellipses to a common collinear motion on the 
degenerate linear ellipse Fb, where the points n and m are 
chosen such that the time of flight (TOF) along nm equals 
the TOF along the elliptical arc NM on the first ellipse. On 
the second ellipse the TOF along the arc NB2M equals the 
TOF along nhm. The points n and m are found by marking 
the point G halfway between A'^ and M. Then the major axis 
Fh = A\B\ — A2B2 of the linear ellipse is drawn starting 
at F and running through G. On this line the point g is 
placed at the distance Fg = ^{FN + FM). Finally n and m 
are given by the intersection points of a circle around g with 
radius GN = GM. This construction shows that the sum of 
the lengths of the shaded triangle a± = FN + FM ± NM is 
equal to a± = Fn + Fm±nm. The ficticious collinear motion 
goes back to Lambert and can be picturized as the limit of an 
elliptic motion with extremely small semiminor a^xis b. The 
eccentricity approaches one from below in such a way that 
the moving particle turns around at F with very high but of 
course still non-relativistic velocity. 



function 5(r,r',T) which follows from the well-known 
Legendre transformation 

S{r,r',T) ^W{r,r',E) - Et . (10) 

The travel time t = t — t' from r' to r can be calculated 
from 

_ dW 

^" ' 

or directly from Eq. ^ by using the method of App. [X] 
The result is 



(11) 




(7 — sin7 — (5 + sin5) . (12) 
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Equation (|12p is Lambert's theorem p. 102) for the 
travel time between N and M. In our case it is more im- 
portant to point out that Lambert's theorem is not only 
valid for the travel time but also for the reduced action 
W and, although not of importance here, for Hamilton's 
principal function S. With these results in mind, we are 
now in a position to calculate the semiclassical Green 
function. 



III. COULOMB TRAJECTORIES AND 
LAMBERT'S THEOREM 




In quantum physics, the Kepler problem becomes the 
Coulomb problem. The connection between classical and 
quantum mechanics is conveniently established through 
the introduction of the quantum mechanical Green func- 
tion, also called propagator. The Green function is the 
mathematical vehicle that allows a particle to go from an 
initial configuration to a final one. In configuration space 
it represents the transition amplitude to travel from r' to 
r. Each classical trajectory in Fig. [T] has sharp energy 
and travel time. In quantum mechanics the travel may 
occur either with fixed energy or in a given time. Travel 
with fixed energy is characterized by the nonrelativistic 
energy Coulomb Green function which was obtained by 
Hostler fl^l in configuration space in closed form, start- 
ing from a partial-wave expansion. The Coulomb Green 
function in momentum space was derived soon after by 
Schwinger [3]. 

Feynman's path integral method is a natural way to cal- 
culate transition amplitudes. For classically allowed tran- 
sitions one has to identify all classically allowed trajecto- 
ries, assign each of them with a phase and an amplitude 
and sum up their contribution. This procedure yields 
as an approximation the semiclassical amplitude. In a 
quantum mechanically exact calculation of the propaga- 
tor one would have to sum up all paths, including the 
classically forbidden ones. Semiclassical methods work 
usually very well because classical trajectories carry the 
main information needed to calculate the transition am- 
plitude from r' to r [l^. In particular, semiclassical 
methods are accurate and useful when large angular mo- 
menta are involved. Typical problems with high-angular 
wave packets require in an exact quantal calculation a 
non-trivial summation over many partial waves of the 
Green function. This problem is avoided in the semiclas- 
sical treatment presented here where no summation over 
partial waves is necessary. For < a particle will move 
on an ellipse in a plane with the center of force in one 
focus. The binding energy fixes the length of the semi- 
major axis while the semiminor axis will also depend on 
the angular momentum. As shown in Fig. [21 the classi- 
cally allowed trajectories are confined to the volume in 
position space defined by the equations of motions for 
a given initial absolute value of the velocity. In two di- 
mensions, one has a critical ellipse that leads to a finite 
classical motion. 




FIG. 2: The caustic (thick solid line) for the Kepler problem 
is an ellipse. Classically allowed orbits with the same energy, 
the same center of force F and a common starting point TV 
have to lie inside this critical ellipse. 



From the definition (1) of a+ and a- it follows that a+/2 
and a_/2 are the distances of m and n from point F (see 
Fig. [1]). Hence we identify the distances a± as path co- 
ordinates of n and m along the straight line Fnm with F 
as origin. Energy conservation H = E in Eq. ^ yields 
the velocities 



v± 



2\E\ /4a -a± 



/i V Oizt 



(13) 



in m and n. By making use of the coordinates a/2 and 
velocities v{a) we obtain the reduced action for travelling 
from n to m: 

W(r,r';E) ^ W+{a+, E) - W-{a-, E) , (14) 

where 

na±/2 

W±=fi v+d{a±/2) 
Jo 



"±/2 Kcyaa±-{a±/2y 



(15) 



can be cast in closed form 



d(ci±/2) 



a V 2 



rrr ^cl^ ( 1 rT, ^ n / CK± 

W± ~ \\ I - V (4a ^ «±)«± + 2aarctan 



4a — a± 
(16) 

The last three equations are consistent with Eqs. ([8]) and 
([9|) and confirm the essentially one-dimensional charac- 
ter of the reduced action for the Kepler and Coulomb 
problem. 



4 



m 




FIG. 3: The four elementary paths from n to m according to 
Lambert's mapping theorem. The travel time for each path 
can be expressed by t_ , t+ and the time for a round trip. 



IV. SEMICLASSICAL ENERGY GREEN 
FUNCTION 



paths i leading from r' to r and having the reduced ac- 
tion Wi. The VVPM-determinant contains derivatives 
of second order with respect to r, r' and E. For ex- 



ample, 



drdr' 



is a n X n matrix with mixed derivatives 



with respect to starting (initial) and ending (final) points 
r' = ix[,X2, ■■■,x'j^) and r = {xi,X2, ■■■,Xn). 
Finally, is the Morse index which is the number 
of conjugate points along the trajectory from r' = 
{x[,x'2, ■■■jx'n) to r = {xi,X2, ■■■,Xn)- In the next sec- 
tion the indices will be read off from the analytical result 
for the VVPM-determinant. We have seen before that 
Lambert's theorem allows the Coulomb problem to be 
mapped on a 1-D problem. Utilizing Lambert's projec- 
tion theorem, we are now in a position to find all possible 
trajectories and, if needed, all traveling times. Figure [3] 
reveals all elementary possibilities to travel from n to to. 
As in Fig. [1] already mentioned, we regard the motion 
in 1-D as motion on an ellipse with infinitesimally small 
semiminor axis b. For such motion we obtain Table [H 
where 



The n-dimensional energy Green function is a solution 
of the inhomogeneous stationary Schrodinger equation 



[E - H] G(") (r, r';E)^ ^("^ (r - r'). 



(17) 



W2^ = 2Tr^/fiaKc and = 2ttJ (21) 



with (5*-"^ being a delta-function point source in n dimen- 
sions. Different boundary conditions on G*-"^ are pos- 
sible. For scattering problems, outgoing-wave bound- 
ary conditions are usually appropriate. For stand- 
ing waves and for bound-state problems G^"^ is real. 
G*^"^ (r, r'; i?) characterizes the probability-amplitude for 
traveling from r to r' with a given energy E. For the n- 
dimensional Coulomb problem the Hamiltonian is given 

by 



denote the action for a closed orbit and the time of circu- 
lation on the same closed orbit respectively. We observe 
that both quantities depend only on the semimajor axis 
a, i.e. on the orbital energy E. 

Traveling from r' to r at constant energy is possible along 
one of the four elementary paths. However there is an in- 
finite number of possibilities for traveling due to the ad- 
dition of an arbitrary number of loops to each elementary 
path. 



H 



2m 



(18) 



where A is the Laplace operator and r the distance from 
the force center in n dimensions. 

The semiclassical limit of the energy Green function is 
given by 



G^-\vy-E)^^Y. 



-1 



ifi ^ (-27rin)("-i)/2 
X exp 



I D(iy,(r,r' ;£;))! 



1/2 



iw^,(r,r';£;)-im,^ 



where 



D(W(r,r';i;)) =det 



drdE 



(19) 



(20) 



dEdr' dE^ 



is the Van Vleck-Pauli-Morette (VVPM) determinant. 
In Eq. (fT9|) one has to sum over all classical fixed-energy 



THE VAN VLECK-PAULI-MORETTE 
DETERMINANT 



We will now calculate the amplitude of the Green's 
function, i.e. the VVPM determinant D(W^) (Eq.dOl)- It 



is helpful to realize that the n x n sub-determinant 
vanishes (j^, page 24). Therefore D(M^) is independent 



of the matrix element 



We will replace this element 



by 0. The VVPM matrix contains mixed second deriva- 
tives of W with respect to the coordinates x'^ and Xj . In 
the last section we showed how to express the reduced 
action for the elementary four paths as combinations of 
the two basic actions W+{a+{r,r')) and Vl^_(Q;_(r,r')) 
(Fig. [3] and Table |l| . By using g(^'^/9) = together 
with the chain rule we find for the off diagonal {i ^ j) 
elements 
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path 



action 



travel time 



Morse index 



® direct path Wi = W+ - W- Ti ^ t+ - 

(2) reflection at f W2 = W+ + W- T2=t++t- 1 

(3) two reflections W3 = W2n + {W+ - W-) T3 ^ T2,r - {t+ - t-) 2 
@ reflection at the caustic W4 = W27, ~ {W+ + W_) T4 = T2,r 1 

TABLE I: Reduced action and Morse- indices rrii for bounded motion along the four elementary trajectories in classically allowed 
regions for an attractive Coulomb potential in three dimensions (see also Fig. [T|. The reduced actions are combinations of W+ 
and W- and the Morse indices in n dimensions are obtained in Sec. [V] 



dxidx': 4(9(q;+/2) \ r' 



(22) 



d^W^ dv- 



dx.dx'^ 4 9(a_/2) \ r' 



■^j -^j -^j 



Xi Xi + X,^ 

r s 



flV- 



{xi Xf){xj x) 



(23) 



where we made use of the fact that according to Eq p3|) the v± arc functions of a± . The diagonal elements follow 
in a similar fashion 



d'^W+ ^ dv+ 



dxjdx'- 4 d{a+/2) \ r' 



3 ~'~ 



2s 



1 - 



(24) 



dxjdx'- 4 9(a_/2) \ r 



r 



(Xj 



(25) 



The mixed derivatives with respect to energy and coor- 
dinates are obtained by utilizing — t. We also have 



dW± 



— t±. Lambert's projection of the Coulomb prob- 

(26) 



lem to a linear, one dimensional motion implies that 
d'^W± dt± dt± 



d{a±/2)dE d{a±/2) d(a±/2) 
We can therefore write 



dxjdE 



dU 



d{a+/2) 



d{a+/2) dxj 



dxjdE 



2r 



Finally, 



di+ d{a+/2) 



dx[dE d{a+/2) dx[ 



92 W_ 



dx'SE 



2r' 



1 



2r 



2s 



2r' 



2s 



1 

v± 



2s 



(27) 
(28) 



2s 



(29) 
(30) 



In principle it is possible to evaluate the second deriva- 
tives for elliptic motion for all points N and M. However 
this is a tedious task. Fortunately Lambert's theorem 
tells us that elliptical Kepler motion can be mapped on a 
degenerate ellipse where motion occurs on a 1-D straight 
line. We use this mapping and assume the coordinate q 
to run along this line from q' = a_/2 to q = a+/2, i.e. 
from point n to point to. In n dimensions we have n — 1 
coordinates X2, X3, . . . Xn that are orthogonal to the tra- 
jectory. Along the trajectory we have Xi = x'^ — for 
i > 2. If we therefore evaluate Eqs. (HHl-dSni) for i,j > 2, 
we observe that the right hand sides vanish except for 
the diagonal matrix elements 



and 



dxjdx'j ' 



"27 



J > 2 , (31) 



J > 2 . (32) 



Each direction i > 2 orthogonal to the straight-line 
trajectory contributes with the same dimensional fac- 
tor J-. Putting everything together we can cast the 
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(n + 1) X (n + 1) VVPM determinant in a simple form: 



D(W) = det 



a(a + /2)a(a_/2) ^ 

J" 



V 





d{a_/2)dE 











d{a+/2)dE 





(33) 

From Table U] we infer that the action W needed for the 
four elementary paths is always a linear combination of 
W+ and The necessary determinants D(W+ ± W^) 
are obtained from Eq. ([33]) with T replaced by ± T- . 
Recalling that W+ is a function («-) only, we 

conclude that 



9(a+/2)9(a_/2) 



0. 



(34) 



Therefore the entry on the top left of VVPM-matrix 
vanishes. The determinant is now easily calculated via 
Laplace expansion. The result is 



T>{W) = 



d{a+/2)dEd{a^l2)dE 



(35) 



We should point out that Eq. ([35)) is also valid for scat- 
tering states if action and velocities are adapted to un- 
bounded motion. 

We now determine the Morse indices mi which are given 
by the order of the zeros of the determinants D(T) along 
path number i. Here we restrict ourselves to the three- 
dimensional Coulomb problem, n = 3. By inspecting 
Fig. [3] we see that on path the velocities are different 
from zero because we have assumed that neither point 
m nor point n is lying on the caustic 6. Hence we have 
mi = 0. On path @ the velocity vanishes at the reflec- 
tion point b. There a pole of first order is generated in 
the determinant. As a result we have m^ = 1, indepen- 
dent of n. Path (2) corresponds to elliptic motion with 
infinitesimally small semiminor axis b with the particle 
(planet or electron) moving around F with infinite ve- 
locity, V oo. Along this path it therefore encounters 
a pole of order n — 2 = 1 at _F, meaning that the parti- 
cle picks up the Morse index m2 = 1. Obviously we have 
7713 = 7712+7714 = 2. Finally by inspecting Eq. (j36p we ob- 
serve that a full round trip picks up an additional phase 
which originates from closing the loop with = w_ and 
s = 0, giving rise to a pole of order ri — 1 = 2 in the de- 
terminant. In other words, closed orbits pick the phase 
m2TT = 2(77 — 1). 



A straightforward evaluation of Eq. ([35]) yields simple 
results for the determinants of the four elementary paths: 



D® = B{W+ - W-) = 



B® = U{W+ +W-) = — 
= -D(3) . 



(n-l) 



(36) 



(n-l) 



VI. E <0: THE BOUND-STATE GREEN 
FUNCTION 



Having found the amplitudes, reduced actions and the 
correct phases we are in a position to evaluate the semi- 
classical Green function in analytic form. We showed 
before that Gi^\r,r';E) consists of the amplitudes for 
the four elementary trajectories plus a summation over 
all possible loops for each elementary path. Therefore we 
can write 



Gi-J (r, r'; E) ^ V\r>{Wdr,r'; E) + j W.^E))] 



^{W,{r,r';E)+jW2AE))~i^{m,+jm2^) 



X exp 

=Gcicm(r,r';-B) x Pgiob(Vp2^, 
I 



(38) 



as a product of the elementary four-path Green func- Green function. The summation over the infinite number 
tion Golem (r, r'; ^^) and and a factor Pgioh{W2Tr , n) = 
^-'^P [y ~ f "^2ii-)] which accounts for the loop 
summation. The factorization is possible because 
T>{Wi{r,r';E) + jW2^{E)) is independent of W2^iE). 
Each loop adds the same non-negative phase to the 
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of loops can be carried out. We obtain 



(a) Quantum mechanical solution 



^exp 
1 i 



27rij 



■ cot 



\2TTh 4 

_ TO27r 
2TTh 4 



(39) 



The poles of Pgiob(W^27r, yield the energy eigenvalues 
of the hydrogen atom. Obviously, they are obtained from 
the poles of the cotangent given by the non-negative in- 
tegers, - ^ = 0, 1,2, .... Using Eq. dH]) together 
with m2ir = 2(n — 1), it is now easy to extract the ex- 
act energy eigenvalues for the hydrogen atom in n > 1 
dimensions 11411 



1 



2h^ (fc + (n-l)/2)2 



[A; = 0,1, 2,, 



(40) 



We should point out that the correct quantization rule 
for the action in n dimensions 

71 — 1 

1^2. = ;i(fc + __) [fc = 0,l,2,...] (41) 

is an integer multiple of h only for odd values of n. 
The elementary four-path Green function Gdcm can be 
written in a more compact fashion because paths lying on 
the same ellipse have the same amplitude as can be seen 
from Fig. [U and Fig. O Their Morse indices are related 
to each other through nij — [n — \) — rrii. Therefore we 
can merge paths ® and (3) and paths (2) and @ pairwise 
together. Then the elementary four-path Green function 
shows a two-path interference pattern. 
Putting everything together, we recast the (real) negative 
energy Green function (|38p in the form 

G(?Hr,r';£;) ^ ^ 



cos 



sm 



i)/2 sin(fc7r) 
Wi 

3(n- 1) 



1 



W2 

h 



(42) 

The bound states (j40|) at k ~ 0,1,2,3,... give rise to 

poles in Gs"'*. Note that in Eq. ((42|l k can assume any 
continuous value fc > 0. DQ and D(2) are the Van Vleck- 
Pauli-Morette determinants given before in Eqs. ((36)) and 
([37ll. The actions Wi = W+ - and W2 = W+ + W- 
are readily calculated from Eq. . Equation is the 
main result of the paper. In the next section we compare 
the semiclassical result for the Green function with the 
exact quantum result. The case E > will be treated in 

App.m 



VII. RESULTS AND COMPARISON WITH THE 
EXACT PROPAGATOR 

In a last step we compare the analytic expressions for 
the Coulomb Green function with the exact quantum 




X (Bohr radii) 

(b) Semiclassical solution 



1000 i 1500 
N 




X (Bohir radii) 

(c) Semiclassical solution with uniform approximation 
1000 




X (Bohr radii) 



1000 i 1500 
N 



FIG. 4; Contour plot of |G*^'(r, r'; E)\'^ with 'principal' quan- 
tum number u — 29.2 and center of force F at the origin. The 
initial position vector N a.t r' = (1232, 0, 0) Bohr radii is lo- 
cated on the X axis, the final position vector r is varied in the 
X, y-plane. 



mechanical Green function as a function of r. We use 
atomic units. Fig. 2] shows contour plots of the three di- 
mensional Green function G{r,r' ; E). In order to avoid 
the infinities at integer principal quantum numbers, we 
have chosen the non-integer 'principal' quantum number 
1/ = k+1 = 29.2 in Eq. (HDl) . This value is close to the one 
treated numerically in [l^ . The center of force is located 
at the origin, the starting point with r' — (1232,0,0) 
Bohr radii has been chosen to lie on the x axis. The end 
point r is varied in the x, y- plane. 

To illustrate the meaning of G we assume to have a co- 
herent stationary source cr(r') of independent particles. 
Such a source will generate the following wave function: 



V'(r, E) 



j G(r, r';E)a{r') dr'. 



(43) 



For a point like source at r' the plot of the Green func- 
tion reveals how particles leak out of the point source at 
r' under the influence of the Coulomb field. In our case 
G is real; hence there is no net current flowing out of the 
source. All particles are eventually reflected back into 
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the source. A comparison with the exact Green func- 
tion shows that all features, including the nodal struc- 
ture are mirrored perfectly by the semiclassical Green 
function. However we must face the fact that the semi- 
classical approach will fail at the caustic where two tra- 
jectories merge into one. Here the deficiency can be re- 
paired by making use of the uniform approximation (see 
App. [C|. The uniform approximation is slightly more 
complex than the semiclassical approximation. 
To demonstrate how well the approximation works we 
present a cut of the Green function parallel to the x-axis 
(Fig. [5]). The semiclassical approximation starts to de- 
viate from the exact solution near the caustic where the 
saddle point approximation that underlies the semiclas- 
sical theory is no longer valid. The spike in the figure 
marks the position of the caustic. At the caustic the 
semiclassical approximation should be replaced by the 
uniform approximation which is seen to match the exact 
quantum solution very well. 

The mapping of the Coulomb problem to a 1-D-problem 
has the great advantage that tunneling properties in a 
Coulomb field can be easily calculated in semiclassical 
approximation because one can avoid the inherent difh- 
culties associated with multidimensional tunneling. Tun- 
neling trajectories are shown in Fig. [71 In the tunneling 
region there is exponential decay but no reflection. The 
analytic continuation of the action into the tunneling sec- 
tor is given in Appendix [Cl The same projection method 
as before can be used. This time the Morse indices are no 
longer integers and will depend on how deep the particle 
will move into the tunnel. 



VIII. CONCLUDING REMARKS 

Lambert's theorem has proven powerful for calculating 
the semiclassical Green function (Eq. ([H])) because it al- 
lows to parameterize all dynamical quantities in terms 
of distances. This feature allowed us to eliminate the 
eccentric anomaly which is ambiguous in the angles and 
therefore has to be treated very carefully The n- 

dimensional Coulomb problem could be reduced to one 
dimensional motion. The reduction is achieved by the in- 
troduction of new variables a± = r + r' ±s. All necessary 
reduced actions could then be found analytically. We 
derived a closed expression for the semiclassical Green 
function. The Morse indices followed directly from the 
analytic form of the Van Vleck-Pauli-Morette determi- 
nant. 

The semiclassical energy Green function is found to be an 
excellent approximation to the exact Green function. It 
also yields the correct bound-state energies for hydrogen 
in all dimensions. We should point out that the semi- 
classical approximation also works very reliably even at 
low energies with small principal quantum numbers. We 
found that also in those cases semiclassics matches the 
quantum mechanical Green function extremely well. For 
small quantum numbers one has less nodes and the ellip- 
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FIG. 5: Plot of G'^' (r, r'; £) with principal quantum num- 
ber u = 29.2 and center of force at the origin. The initial 
position vector is located at r' = (1232,0,0) Bohr radii, the 
final position vector shown along the x-axis for fixed y = 400 
Bohr radii, (a) Quantum mechanical Green function, (b) de- 
viation of the semiclassical Green function from the quantum 
mechanical one, (c) deviation of the uniformly approximated 
Green function from the quantum mechanical result. The 
uniform approximation agrees with the quantum mechanical 
result better than 1/100. 



tically shaped caustic shrinks. 

In energetically forbidden regions there are no classical 
trajectories. Nevertheless we can continue the semiclassi- 
cal Green function into the tunnel. The exit of the tunnel 
can be dealt with in semiclassics by invoking corrections 
given by the uniform approximation. 
The motion in a Coulomb potential is an important prob- 
lem in quantum mechanics. It is therefore useful to learn 
how the semiclassical limit of the energy Coulomb Green 
function emerges from a coherent summation of all am- 
plitudes that belong to an infinite number of classical 
trajectories. The results of this paper can be readily im- 
plemented into real-space problems in the presence of 
Coulomb interaction. One example is the quantum be- 
havior of Rydberg atoms [l^, [l^ . 
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APPENDIX A: REDUCED COULOMB ACTION 

In order to eliminate the eccentricity e in Eq. (|14p in 
favor of spatial positions we introduce in a first step the 
new variables 

cos(7:=ecosl — - — I and a := — - — . (Alj 

In a second step we substitute 

J := d + g and S := g - d (A2) 

to arrive at Eq. ^ . We next relate the variables 7 and 6 
to spatial positions. The radial position of any point M 
or N on the ellipse (see Fig. [T]) relative to the center of 
force is given by 

r^ = [{x-eaf+y% (A3) 

where x and y are the coordinates relative to the center 
of the ellipse. With the help of Eq. ^ we easily find 

r = a(l-ecos^) (A4) 



with Wi^W+- 



2. Repulsive Coulomb Interaction 

If the potential is repelling we have again two hyper- 
bolic trajectories which connect TV and M. But now a 



In terms of the variables g and d we have 

r + r' = 2a(l — cosgcosd) (A5) 

and 

|r — r'l = 2a| sindcos^l (A6) 

Without loss of generality we can assume < d < tt 
and < g < Tr/2. From the last two equations we then 
readily confirm the desired result, Eq. ([9]). 



APPENDIX B: E>0: SCATTERING STATES 

To treat scattering states in semiclassical approxima- 
tion we can use the same formalism as for bound states. 
In an attractive force field and for E > 0, there is no 
caustic and hence no reflection at b. As can be seen 
from Fig. [6] we then have only two hyperbolic trajecto- 
ries leading from N to M. The one-dimensional variables 
are again a± = r + r' ± |r — r'|. The projection of the 
motion onto a line applies again but we have to consider 
the change in geometry. 



1. Attractive Coulomb Interaction 

In this case we obtain 



(Bl) 



(B2) 



I 

caustic separates the classically allowed region from the 
energetically forbidden region. Classically allowed mo- 
tion occurs for 4|a| < a±. The corresponding reduced 
action reads 



a 



(n) 



,(r, r'; £;) = -- 



n (27ri;i)("-i)/2 



V^D^exp[^M/i] + 



\eM^W,~i^{n~2)]] , 



W- and W2 = W+ + W-. The action follows again from Eq. p^ . adapted to £^ > 0, 



(M^±)attr = 




MO J2a 



«±/2 



/2a + a±/2 
a±/2 



d(5±/2) 



— -VA4a+^^^)a^+2alog^^ )= 

a \Z y/Aa J 
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FIG. 6: For E > 0, hyperbolic trajectories connect A'^ and M, with the center of force at F. In the repulsive case (not shown) 
there is a caustic in contrast to the attractive case where every point in space can be reached. 



5±/2 



d(fi±/2) 



iKcfi ( 1 / ^ Vol - v/-4|a| +a± 



For completeness we write down the action in the classically forbidden tunneling region where a_ < 4|a| 



(B3) 



/ 2|a| -5-/2 
V 5-/2 



d(<5_/2) 



= ±i^/^P ^-7r|a| + iV(4|a| - a_)a_ + 2|a| arctan 



(B4) 



4|a| - a- 



APPENDIX C: ANALYTIC CONTINUATION 
INTO THE TUNNELING REGION 

1. Uniform Approximation 

The semiclassical Green function is derived from the 
exact expression for the quantum mechanical Green func- 
tion by making use of the saddle point approximation 
(SPA). However this approximation is not valid at the 
caustic where two saddle points merge into one. In this 
case the uniform approximation (UA) will cure the de- 
ficiency of the SPA. The method is standard. For more 
details the reader is referred to 'F], p. 118ff, p. 131ff. Here 
we follow the method outlined in reference [TTj] . For n = 3 
we have 



and 



2/3 



if %W- 
if ^{W- 



= 0, 



(C3) 



>l/4 



v/D(VK+) + ^-T>{W- 



3— 5i7r/4 



-1/4 



^— 5i7r/4 



(C4) 



2. Tunneling Regime 



Gua(r,r';i;) 



with 



(doAi(-C) - i'^iAi'(-C)) (CI) 



.{lKiW+ + W.) if3(W^_) = 0, 

\^{W+/fi) if 3(W^-) = -3(W+), 



Here we look at classically forbidden motion > 4a. 
This means that point m is lying in the tunnel (see 
Fig. [7]). Whereas W- is not changed, W+ must be contin- 
ued into the classically forbidden space sector. The ana- 
lytic continuation is obtained by making use of Eq. ()16p 



11 



{W+) 



forb 



'A' 

^<■\l — 



2a 



«±/2 



1 



'2a + a±/2 
^/4^ 



d(ci±/2) 



(C5) 



■±i oV(-4a + a±)a± +2alog^= ^ — ^ 



r 




FIG. 7: Bound states and tunneling trajectories: We use the 
projection formalism to find two possible paths to a point 
in the classically forbidden region. As the electron passes 
through b into the tunneling space sector, action and velocity 
become complex (see Eq. (ICSp ). 

There are two complex conjugated solutions. For the 
propagator we select the term with positive imaginary 
part to ensure that the wave function decays exponen- 



tially deep in the tunnel. Note however that both solu- 
tions will contribute to the uniform approximation in the 
vicinity of the tunnel exit. 
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